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ABSTRACT 

Knowledge of late K and M dwarf metallicities can be used to guide planet searches and constrain 
planet formation models. However, the determination of metallicities of late-type stars is difficult 
because visible wavelength spectra of their cool atmospheres contain many overlapping absorption 
lines, preventing the measurement of equivalent widths. We present new methods, and improved 
calibrations of existing methods, to determine metallicities of late-K and M dwarfs from moderate 
resolution (1300 < R < 2000) visible and infrared spectra. We select a sample of 112 wide binary 
systems that contain a late-type companion to a solar-type primary star. Our sample includes 62 
primary stars with previously published metallicities, as well as 50 stars with metallicities determined 
from our own observations. We use our sample to empirically determine which features in the spectrum 
of the companion are best correlated with the metallicity of the primary. We find ~ 120 features in K 
and M dwarf spectra that are useful for predicting metallicity. We derive metallicity calibrations for 
different wavelength ranges, and show that it is possible to get metallicities reliable to < 0.10 dex using 
either visible, J-, H-, or K- band spectra. We find that the most accurate metallicities derived from 
visible spectra requires the use of different calibrations for early-type (K5.5-M2) and late-type (M2- 
M6) dwarfs. Our calibrations are applicable to dwarfs with metallicities of —1.04 <[Fe/H]< +0.56 and 
spectral types from K7 to M5. Lastly, we use our sample of wide binaries to test and refine existing 
calibrations to determine M dwarf metallicities. We find that the £ parameter, which measures the 
ratio of TiO can CaH bands, is correlated with [Fc/H] for super-solar metallicities, and £ does not 
always correctly identify metal-poor M dwarfs. We also find that existing calibrations in the K and 
H bands are quite reliable for stars with [Fe/H]> —0.5, but are less useful for more metal-poor stars. 
Subject headings: binaries: visual - stars: abundances, - , stars: fundamental parameters, - stars: 
late-type, - techniques: spectroscopic 



1. INTRODUCTION 

Despite their intrinsic faintness, M dwarfs ha ve become 
attractive targets for exoplanet searches 



neau et al 
M dwarts 



2009 Apps et al. 2010 Mann et aT'20TT), 



Charbon- 



have smaller radii and lower masses than 



their solar-type counterparts, allowing for easier detec- 
tion of low- mass exoplanets (Gaidos et al. [2007). For 
transiting M dwarf planets with hig h-precision observa- 



tions, such as those d one by Kepler ( Borucki et al. 112010 



Batalha et al. 2012), errors in planet parameters (pri 
marily planet rad ius) are directly tied to errors in stellar 



parameters (e.g., Muirhead et al. 2012), which in turn 



depend on reliable measurements of stell ar temperature 
and metallicity (potter et al.||2008"| IDemory et al.||2009| 
Allard et al.||20lT ). Accurate metallicities are also nec- 
essary to study any c orrelation between metallicity and 
plane t frequency (e.g. Johnson et al. 2010 Mann et al. 
2012|. Further, because the stellar mass function peaks 



around mid-M, M dwarfs weigh heavily in any study of 
Galactic structure. The distribution of M dwarf metal- 
licities can therefore be used to set limits on Milky Way 
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formation models ( Woolf fc West|2012 ). However, unlike 
those of their solar-type counterparts, M-dwarf metal- 
licities are difficult to determine, primarily due to the 
presence of complex molecular lines in their visible spec- 
tra, which result in line confusion and a lack of iden- 
tifiable continuum, an d do not always match with cur 



rent M dwarf models (Mould 
The visible wavelengths usual 



1976 



Allard et al. 2011). 



y used to derive metallic!- 



ties for solar-type stars are dominated by TiO lines that 
increase in strength with decreasing effective tempera- 
ture. These TiO lines obscure the continuum, making 
equivalent width measurements unreliable, and differ as 
a function of spectral type. Direct spe ctral synthesis has 



been tried on a sma ll sample of stars ( Bean et al. 2006 
Onehag et al. 2012). However, such spectra is observa- 
tionally expensive, and spectral synthesis is complicated 
by incomplete lines lists. 

One common technique is to use wide binaries with 
a solar-type (late-F, G, or early-K) primary and a late- 
type (late K or M) dwarf companion. Since these stars 
presumably formed from the same molecular cloud, the 
metallicity of the M dwarf companion can be a ssumed 
to be the s ame as that of the FGK primary (Bonfils 
et al. 2005). From this one can create empirical cali- 
brations of observable features to determine metallicities 
in M dwarfs. Techniques based on absolute photometry 
dBonfils et al.|2 005|| Johnson fc Apps|2009[ iSchlaufman fc 
Laughiin||2010| 



aughlin 2010) require parallaxes, which are only avail 
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able for a handful of the closest and brightest M dwarfs 



Molecular indices at visible waveleng ths (e.g. Woolf fc 



Wallerstein 2006 Lepine et al. 2012[ ) have proven use 
tul tor separating M dwarfs into luminosity / metallicity 
classes (i.e., dwarf, subdwarf, and extreme subdwarf), 
but saturate near solar metallicity and are less reliable 
for late-K and early-M dwarfs. 

The use of spectral indices in the infrar ed has, how- 
ever, been s howing considerable promise. |Rojas-Ayala| 
et al. (2010 henceforth RIO) showed that [Fe/HJ can be 
interred tor M dwarfs using the Na I doublet and Ca I 
triplet in moderate resolution (R ~2700) i f -band spec- 
tra and assu ming solar relative abundances. |Rojas-Ayala| 
et al. (2012 henceforth R12)) expanded on this, showing 
that one can determine [M/H] for early to mid-M dwarfs 
as ac curately as 0. 1 dex using the same indices. Fur- 
ther, Terrien et al. ( 2012| henceforth T12) were able to 
demonstrate similar precision in the H— band. However, 
calibrations for both H- and if-band metallicities were 
derived using a relatively small sample (18 and 22, re- 
spectively) of wide binaries. As a result, they were only 
verified for systems with near-solar metallicities (—0.4 < 
[Fe/H] < +0.3), and for a narrow range of spectral types 
(M0-M4). 

In this paper we analyze visible and near-infrared 
(NIR) observations of 112 late-K and M dwarfs with 
F, G, or early K star primaries. We provide determine 
metallicity estimates for 50 of the FGK stars. We also 
provide a full list of the 120 features we identify to be 
metal-sensitive features in visible and NIR spectra, with 
the expectation that these will be useful for future stud- 
ies on M dwarf metallicities. In Section [2] we present our 
sample of wide-binary systems, followed by a descrip- 
tion of our observations and reduction in Section [3l In 
Section [4] we explain how we determine basic properties 
for both the primary and companion stars. We discuss 
our search for metal-sensitive features in Section [5] and 
present the results from this analysis in Section [6j which 
includes our calibrations to determine metallicities at a 
range of wavelengths. We test existing techniques in Sec- 
tion [7J Lastly, we summarize the results in Section 8 
and discuss possible drawbacks of our analysis as well as 
prospects for future studies. All wavelengths used in this 
work are stated as vacuum values. 



2. SAMPLE 

We draw a sample of late K and M dwarfs with co 

moving F GK stars from a variety o f literature sources, 

s pecifi ca lly|Chaname fc Gou ld (2004),|Gould & Chaname 

£1 d'vi/r 



specmcairy onanamc & (jould (z!UU4), uo 
P004l), iLejise fc Bo^gioTnop[2W|~ 
Laughlin (20lU|), and Tokovinin & Lepi 
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Schlaufman & 
We" 



(2012). 



also identity a number ot new common-proper motion 
(CPM) pairs that contain an FGK star and a late K or M 
star using proper motions from (in order of preference 



star using proper motions trom ( m order ot prctcrence ) 
Hippa rcos (|van Leeuwen fc~Fantino 2005 van Leeuwen 



SUPERBLLNK (Lepine & Shara 2005), or the 



XL survey ( |Roeser et al.||20r0| r We identify new 
CPM pairs, as well as vet the literature sample, f ollow- 
ing the techniques of Lepine fc Bongiornol (|2007[ ) and 
Dhital et al.] ( |2010| ). |Lepine fc Bongiorno| ( |2007f ldentify 
CPM pairs using the formula: 



where /i is the absolute proper motion of the primary star 
in arcsec yr -1 , AO is the separation between the two stars 
in arcsec, and A/i is the absolute difference in proper 



Lepine 



motion between the two stars in arcsec yr -1 
& Bongiorno (2007) consider two stars to be physically 
associated with each other if AX < 1, AO < 1500", and 
A/i < 100 mas yr -1 . We take a more conservative cut 
and require that our pairs have AX < 0.9, AO < 650", 
and A/i < 60 mas yr -1 to cut down on the number of 
chance alignments. This technique takes advantage of 
the absolute proper motion, and not just the difference 
in proper motion, of the star, however, it docs not factor 
in associated errors in p roper motion. Thus we further 
force the constraint from Dhital et al. (20101: 



A/iq 



( ^Hs 



< 2, 



(2) 



where A/i a and A(j,$ are the differences in proper motion 
between the two components (right ascension and decli- 
nation, respectively) and <ta/j is the error in the proper 
motion differences. Equation[2]is designed to ensure that 
the resulting sample has minimal contamination from 
chance alignment pairs, but at the cost of excluding a 
large number of true wide binaries. To increase our sam- 
ple, we add pairs that do not satisfy Equation 2 but have 
published parallaxes for both the primary and companion 
star consistent to la or the parallax uncertainty, which 
strongly suggests that the pair forms a physical system 
We remove pairs with 5 < —45° or S > +68°, since 
these are outside the reach of telescopes in Hawaii, and 
pairs with AO < 3" as these will be difficult to observe 
and may have inaccurate photometry. We further re- 
move pairs with V c — J c < 2.5 (the subscript c denoting 
the companion and p the primary). This cut will re- 



move almost all stars earlier than K5 (Lepine & Gaidos 
20111 where metallicities can be measured using mod- 
ified spectral synthesis techn iques (Valenti & Piskunov 
1996 Valenti & Fischer 2005 ). We cut out systems with 
V c > 18, and K c > 12, as the observation time required 
for these stars is highly prohibitive. For the same reason, 
we remove systems with V p > 12 unless the primary star 
already has a published metallicity. The resulting sam- 
ple contains 262 pairs with primaries with colors in the 
range 0.8 < V p — J p < 2.5 and/or published te mperatures 
consistent with a late-F , G, or early-K dwarf (Fitzgerald 
1970 Ducati et al. 2001 ). From here we prioritize our ob- 
servations based on; (1) V c and K c magnitudes (to mini- 
mize required telescope time), (2) V c — J c color (to ensure 
a range of spectral types) , (3) availability of metallicities 
for the primary star in the literature or V p magn itud e 
if no literature metallicity is available (see Section 3.1), 
(4) availability of parallax information for the primary 
or secondary, which makes it more likely that these are 
true CPM pairs, and (5) metallicity of the primary (if 
available) to guarantee a range of metallicities for our 
analysis. 

The resulting sample includes 112 stars for which 



AX = [(/i/0.15)- J ' 8 A6A/i] 



10.5 



(1) 



we have obtained infrared spec tra (Section 3.2), visible 
wavelength spectra (Section |3.3| , and metallicities for the 
primary star from the literature or from our own analysis 
(Section [O). Our sample has -1.04 <[Fe/H] p < 0.56, 
and companion star spectral types from K5.5 to M6. We 
list our sample in Table [l] 
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Table 1 



Companion Primary 
Common Name R.A S SpT a Name [Fe/H] 6 ± Source c [M/H] 6 ± Source c 



NLTT 738 


3.8107 


+53.0794 


M2.2 


HIP 


1224 


+0.07 


0.03 


This work 


+0.03 


0.03 


This work 


NLTT 923 


4.6075 


+44.0272 


M4.1 


HIP 


1475 


-0.05 


0.03 


VF05 


-0.53 


0.03 


VF05 


NLTT 2478 


11.3066 


+0.2642 


M3.8 


HIP 


3540 


+0.02 


0.03 


VF05 


-0.01 


0.03 


VF05 


NLTT 3598 


16.3740 


+15.3885 


K7.4 


HIP 


5110 


-0.09 


0.03 


This work 


-0.13 


0.03 


This work 


NLTT 3725 


16.9105 


+22.9558 


M3.9 


HIP 


5286 


+0.25 


0.03 


This work 


+0.21 


0.05 


This work 


Gl 56.3 


19.6676 


—0.8744 


K7.6 


HIP 


6130 


—0.01 


0.06 


C01 


+0.01 


0.06 


Cll 


NLTT 4568 


20.6531 


+12.7524 


K7.8 


HIP 


6431 


+0.10 


0.03 


This work 


+0.04 


0.03 


This work 


NLTT 4599 


20.7525 


-12.9583 


K7.9 


HIP 


6456 


+0.45 


0.03 


This work 


+0.36 


0.06 


This work 


Gl 81.1 


29.2961 


-10.2481 


M0.1 


HIP 


9094 


+0.12 


0.03 


VF05 


+0.05 


0.03 


VF05 


Gl 100 


37.1329 


-20.0407 


M2.9 


HIP 


11565 


-0.28 


0.04 


N12 








NLTT 8107 


37.2792 


+22.8671 


K5.9 


HIP 


11572 


-0.08 


0.03 


This work 


+0.00 


0.03 


This work 


Gl 105 


39.0636 


+6.8717 


M3.9 


HIP 


12114 


-0.12 


0.03 


VF05 


-0.00 


0.03 


VF05 


NLTT 8787 


41.0427 


+49.2317 


M1.6 


HIP 


12777 


+0.06 


0.03 


VF05 


+0.02 


0.03 


VF05 


Gl 118.2 


43.8991 


+26.8724 


M3.8 


HIP 


13642 


+0.28 


0.03 


VF05 


+0.21 


0.03 


VF05 


NN 3195 


46.1810 


+61.7358 


M2.9 


HIP 


14286 


-0.28 


0.03 


Sll 


-0.26 


0.06 


Cll 


NLTT 10349 


48.7538 


+1.0523 


M1.0 


HIP 


15126 


-0.92 


0.03 


This work 


-0.66 


0.03 


This work 


NLTT 11125 


53.0251 


+43.6669 


K7.3 


HIP 


16467 


-0.01 


0.03 


This work 


-0.01 


0.03 


This work 


NLTT 11176 


53.3085 


+46.2553 


K7.3 


HIP 


16563 


+0.20 


0.03 


This work 


+0.16 


0.06 


This work 


NLTT 11270 


53.8689 


+42.8930 


M0.4 


NLTT 11280 


-0.40 


0.08 


Cll 


-0.30 


0.06 


Cll 


NLTT 11500 


54.8955 


+18.3143 


M1.8 


HIP 17076 


-0.52 


0.08 


Cll 


-0.44 


0.06 


Cll 


PM J0355+5214 


58.9037 


+52.2414 


M2.7 


HIP 18366 


-0.36 


0.05 


Ra07 


-0.33 


0.06 


Cll 



Note. — Table^is published in its entirety in the electronic edition of the Astronomical Journal, and can be downloaded 
with the arXiv version of the manuscript. A portion is shown here for guidance regarding its form and content. 

Spectral types derived from TiO and CaH indices, (see |Lepine et ali] ( [20 12 ) ) . Continuous spectral types (to 0.1) are 
cd for plotting/binning/calculations, even though spectral types are only accurate to ±0.2 (and by convention should be 



used for plotting/binning/calculatii 
rounded to the nearest 0.5). 

b [M/H] and [Fe/H] values s hown here include our applied correcti ons (see Section [4 | . 
Metall icity sources: C01 =|Cayrel de Strobel et al.| |200H, M0 4 =|Mis henina et afTl|2004 |l 

l'akeda"et al.l 



VF05 
Ro07 



Valenti &c Fischerl 
tobinson et aLllCTCf 



7t~FD8 = 



05 



B06 = jrjean etalJ 
Ida Silva et al.M^UlI 



LH05 
Ra07 



Luck & Heiter 



TTTT 



(2005), 



Kamirez et af 

|Fuhr mann (200Sjr| !sll =ida Silva et al. (2011), Cll' = |Casagrande et aLl^UT'l'y 
|Neves~et al.| l [2U12| l, I'ms work = analysis ol (Jt til'/ ESPaDOnS spectra as part of this program. 
Note that all metalficities sources are from high-resolution spectra, with the exception of Ro07 
resolution spectra, and Cll, which uses Stromgren photometry. 



(2007), 
^JT2 



vhich uses moderate- 



3. OBSERVATIONS AND REDUCTION 

3.1. ESPaD OnS/CFHT 

Between 2011 January and 2012 April, 60 F-, G- and 
early K-type stars were observed using the ESPaDOnS 
spectrograph attac hed to the C anada France Hawaii 
Telescope (CFHT; 



Donati 20031 on Mauna Kea. Ob- 



servations were taken in the star+sky mode, which gave 
a resolution of R~ 65000 and a wavelength range from 
0.37/im to 1.05/im. All observations were designed to 
achieve a signal-to-noise ratio (S/N) of > 100 at 0.67/im, 
and typical S/N was > 150 (per resolving element). The 
data were reduced au tomatically us i ng the Libre-ESpRIT 
pipeline described in Donati et al.| ( |1997[ ) . 



> 150 (per resolving element). To correct for telluric 
lines, we observed an AOV-type star within 30 minutes 
and 0.1 air mass of the target observation (and usually 
much closer in time and air mass). Often the same A0V 
star was used to remove telluric lines for more than one 
target. To remove effects from large telescope slews, we 
obtained flat-field and argon lamp calibration sequences 
after each A0V star. 
Spectra were extracted and reduc ed using the SpeX- 



Tool package (Cushing et al. 2004), which performed 



3.2. SpeX/IRTF 

We obtained near-infrared spectra of our sample of 
comp anions using the SpeX spectrograph ( Rayner et al. 
|2003b attached to the NASA Infrared Telescope Facility 
(TRTF) on Mauna Kea. SpeX observations were taken in 
the short cross-dispersed (SXD) mode using the 0.3x15" 
slit, yielding simultaneous coverage from 0.8 to 2.4/im 
and a resolution of R ~ 2000. The star was placed at 
two positions along the slit (A and B). Exposures were 
taken with an ABBA slit-nodding pattern, with at least 6 
exposures in total. Integration times were no longer than 
120 s for each exposure to minimize the effect of changes 
in atmospheric H2O. Thus for faint stars more than 6 
exposures were required to get sufficient S/N. S/N in the 
H- and if -bands was > 100 in all cases, and typically 



flat-field correction, wavelength calibration, sky subtrac 
tion, and extraction of the one-dimentional (ID) spec- 
trum. Multiple exposures were stacked using the IDL 
routine xcombxpec. A telluric correction spectrum was 
cons tructed from each A0V star using the xtellcor pack- 
age ( Vacca et al.||2003 ), and then applied to the relevant 
target spectra. 

Reduced spectra we re put in vacuu m wavelengths us- 
ing the formula from Ciddor (1996). We put spectra 
in the star's rest frame by comparing them to a spec- 
trum of the template star HD36395 (an Ml. 5 dwarf, also 
in rest frame/vacuum) taken from the IRTF spectral li- 



brary ( Cushing et al.|2005| |Rayner et al.|2009[ ) . We cross- 
correlated each spectrum with the template, in orders 3 - 
7 separately (order 8 is ignored because it is too smooth 
and has relatively poor S/N), yielding 6 radial velocities 
(RVs). We shifted the spectrum by the average (after 
removing any 5cr outliers) of each set of RVs. 
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3.3. SNIFS/UH2.2m 

We obtained a visible spectrum of each companion 
with the SuperNo va Integral Field Spectrograph (SNIFS, 
Lantz et al.|2004 ) on the University of Hawaii 2.2m tele- 
scope atop Mauna Kea. SNIFS has R ~ 1300 and splits 
the signal with a dichroic mirror into blue (0.32-0.52 fim) 
and red (0.52-0.95 /xm) channels. SNIFS data processing 
is perfor med with a data red uctio n pipeline, described in 
detail in |Bacon et"al] ( |2001[ ) and [Aldering et al.| ( 12006[ ). 
SNIFS processing includes dark, bias, and Hat-held cor- 
rections, assembling the data into red and blue three- 
dimcntional data cubes, and cleaning them for cosmic 
rays and bad pixels. Wavelengths are calibrated with 
arc lamp exposures taken at the same telescope point- 
ing as the science data. The calibrated spectrum is then 
sky-subtracted, and a 1-D spectrum is extracted using 
a point-spread function model. Corrections are applied 
to the spectrum for instrument response, and for telluric 
lines based on observations of the Feige 66, Feige 110, 
BD+2 84211, or BD+174708 spectrophotometric stan- 
dards ( |Oke||199"0 ) that are taken over the course of each 
night. 

Approximate spectral type s are determined by the 
HAMMER software package ( [Covey et aT1 l 2007| ). The 
spectra are then shifted to zero radi al velocity by cross- 
correlating with templates from Bochanski et al. (20071 
of the corresponding spectral type. Late-K stars are 
cross-correlated using an M0 template. 



3.4. Construction of a combined M dwarf Spectrum 

We use the overlapping region in our SpeX and SNIFS 
data (0.81-0.96/xm) to combine the visible and NIR spec- 
tra. We normalize each SNIFS spectrum by a constant, 
C, which is equal to the ratio of the median flux of the 
SNIFS spectra in overlapping region to the median flux of 
the SpeX spectra in the overlapping region. We also find 
that there is a systematic offset in wavelength between 
the visible wavelength and NIR spectra in the overlap- 
ping region. This amounts to a RV shift of ~ 30 km s _1 
between the SNIFS and SpeX spectra. The most likely 
explanation for this is a small difference in the RV tem- 
plates used for our Spe X and SNIFS data, tak en from 
Rayner et al. ( 2009 ) and Bochanski et al. ( 2007 ) respec- 
tively. We choose to shift the visible wavelength data to 
match with the NIR data (which has higher resolution 
and therefore gives more reliable RVs) by adding an addi- 
tional RV correction of —30 km s _ l to each SNIFS spec- 
trum. Given the modest resolution of SNIFS data, this 
correction is unlikely to significantly change our results. 
The offset corresponds to < 1 A, which is significantly 
less than the resolving power of SNIFS, and is similar in 
size to random errors in our RV measurements. 



4. DERIVING STELLAR PARAMETERS 

4.1. FGK Metallicities 

We draw primary star metallicities from a variety of 
sources in the literature. We list the adopted metallicity 
and literature source for each binary in Table [II In total, 
33 p rimary star metallicities come from the SPOCS cata- 
log ( Valenti fe Fischer|2005 ), 50 from observations taken 
as part of this project with CFHT/ESPaDOnS, and 29 
from other literature sources. SPOCS consists of high- 
resolution echelle spectra of > 1000 F-, G-, and K-type 



stars obtai ned with the Keck, Lick , or Anglo- Australian 
Telescope. Valenti & Fischer ( 2005 ) fit the observed spec- 
trum to a synthetic spectrum usin g the software pack- 
age S ME (Spectroscopy Made Easy, Valenti & Piskunov 
1996), which provides a set of observational parameters 
(1'eff , [Fe/H], [M/H], log g, etc.) for each star. We adopt 
an uncertainty of 0.03 dex for their derived [M/H] and 
[Fe/H] values. 

To determine the stellar parameters of primaries ob- 
served with CFHT/ESPa DOnS we model each spec trum 
using the SME software ( Valenti fe Piskun"ovj|1996[ ), fit- 
ting the spectrum to a set of tuned lines from the SPOCS 
catalog (Valenti & Fischer 2005). We simultaneously 
solve for surface gravity, effective temperature, projected 
rotational velocity, and metallicity in addition to individ- 
ual abundances of Fe, Na, Si, Ni, and Ti. Solar values 
are assumed for all of the initial models and after ob- 
taining an initial fit, we then perturb T e g- by ±100K and 
fit again. Our final model parameters are x 2 -weighted 
averages of three runs. Corrections base d on Vesta and 
stellar bina ry observations as detailed in [Valenti fe Fis-| 
cherl(p005| are then applied. 



For stars with good parallax measurements (Hippar- 
cos stars or their companions), w e use Yonsei-Yale (Y 2 ) 
isochrones ( jDemarque et al. | 2004h t o better constrain the 
surface gravity ( | Valenti et al.||2009| . After we determine 
the stellar parameters as above, we use distance and B 
and V magnitudes to the derive bolometric luminosity 
and, combined with the SME T c ff , the stellar radius. 
Bolometric corrections are obtained by interpolating in 
the high temperature grid of VandenBerg & Clem ( 2003), 
and B, V magnitudes we re drawn from hipparcos 



van 



Leeuwen & Fantino 2005 1 . The SME determined ratio 
of Si to Fe is used as a proxy for alpha element enhance- 
ment. A best-fit evolutionary model is found by inter- 
polating in the Y 2 grid which yielded a surface gravity 
for the star. This log g is compared to the value deter- 
mined using SME and if the two did not match, a new 
set of SME models is found with the gravity fixed to the 
isochrone value. The process is repeated until the log g 
values agree to within 0.001 dex. Final stellar parame- 
ters (T e g , log g, [M/H], etc.) for stars observed as part 
of our program are listed in Table [2] 

Our analysis of the ESPaDOnS spectra is designed 
to keep our metallicities consistent with those from the 
SPOCS catalog (both are based on SME analysis and use 
the same set of spectral lines) . As an extra check on con- 
sistency we have obtained CFHT spectra for three stars 
in the SPOCS sample. The derived stellar parameters 
from these three spectra are consistent (within errors) 
with those listed in the SPOCS catalog, confirming that 
there is no systematic offset between metallicities from 
SPOCS and CFHT. 

Metallicities from other literature sources are not nec- 
essarily determined in the same way as our spectral anal- 
ysis, and thus may have small systematic inconsistencies. 
We correct for this by checking for overlap between the 
SPOCS samples and any given literature source. For 
us to use a metallicity derived from any other literature 
source we require; (1) at least 30 stars in both sam- 
ples that can be used as a control sample to check for 
differences, (2) metallicities for our primary stars from 
the literature source fall within the range of metallicities 
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Table 2 

Parameters of Primary Stars observed at CFHT 



Name 



T cff ± log g ± [Fe/H] ± [M/H] ± [Na/H] ± [Ti/H] ± [Si/H] ± [Ni/H] ± 



~y 1 ~ 

A rnd 



Run Type" 



HIP 1224 5141 44 4.53 0.06 +0.07 0.03 +0.03 0.03 +0.12 0.03 +0.03 0.05 +0.02 0.02 +0.08 0.03 2.9 ITER 

HIP 5110 4648 44 4.61 0.06 -0.09 0.03 -0.13 0.03 -0.25 0.03 -0.18 0.05 -0.15 0.02 -0.11 0.03 6.0 ITER 

HIP 5286 4676 53 4.60 0.06 +0.25 0.03 +0.21 0.05 +0.30 0.09 +0.14 0.05 +0.09 0.06 +0.26 0.03 14.7 ITER 

HIP 6431 4858 44 4.57 0.06 +0.10 0.03 +0.04 0.03 +0.10 0.03 +0.00 0.05 +0.06 0.02 +0.11 0.03 6.3 ITER 

HIP 6456 5222 52 4.44 0.06 +0.45 0.03 +0.36 0.06 +0.62 0.05 +0.27 0.08 +0.43 0.02 +0.53 0.03 6.8 ITER 

HIP 11572 5093 44 4.51 0.06 -0.08 0.03 +0.00 0.03 +0.05 0.03 +0.03 0.05 +0.05 0.02 -0.02 0.03 3.6 ITER 

HIP 15126 5285 44 4.64 0.06 -0.92 0.03 -0.66 0.03 -0.84 0.03 -0.56 0.05 -0.51 0.02 -0.86 0.03 2.2 ITER 

HIP 16467 5539 44 4.33 0.06 -0.01 0.03 -0.01 0.03 -0.09 0.03 -0.04 0.05 +0.02 0.02 -0.04 0.03 2.1 ITER 

HIP 16563 5788 44 4.52 0.06 +0.20 0.03 +0.16 0.06 +0.10 0.03 +0.21 0.05 +0.04 0.02 +0.05 0.04 4.4 ITER 

HIP 31127 5158 44 3.83 0.06 -0.54 0.03 -0.45 0.03 -0.56 0.03 -0.50 0.05 -0.34 0.02 -0.56 0.03 3.1 ITER 

HIP 31597 5428 44 4.44 0.06 +0.09 0.03 +0.08 0.03 +0.14 0.03 +0.08 0.05 +0.10 0.02 +0.12 0.03 2.0 ITER 

HIP 32423 4817 44 4.64 0.06 -0.26 0.03 -0.21 0.04 -0.32 0.04 -0.27 0.07 -0.24 0.03 -0.27 0.03 5.6 ITER 

HIP 35449 6156 44 4.36 0.06 +0.21 0.03 +0.18 0.03 +0.17 0.03 +0.21 0.05 +0.19 0.02 +0.18 0.03 1.9 ITER 

HIP 40298 5618 44 4.51 0.06 -0.07 0.03 -0.09 0.03 -0.09 0.03 -0.04 0.05 -0.08 0.02 -0.14 0.03 1.7 ITER 

NLTT 12373 5838 93 4.62 0.09 -0.08 0.05 -0.00 0.08 -0.08 0.05 +0.05 0.08 -0.04 0.03 -0.12 0.06 2.2 VESTA 

HIP 45863 5172 44 4.51 0.06 -0.12 0.03 -0.13 0.03 -0.07 0.03 -0.10 0.05 -0.09 0.02 -0.12 0.03 4.3 ITER 

NLTT 23002 5259 44 4.58 0.06 -0.15 0.03 -0.11 0.03 -0.12 0.03 -0.11 0.05 -0.14 0.02 -0.16 0.03 3.1 VESTA 

HIP 50802 4472 44 4.70 0.06 -0.01 0.03 -0.06 0.03 -0.25 0.04 -0.11 0.05 -0.21 0.06 -0.14 0.03 13.1 ITER 

HIP 54155 5547 50 4.54 0.06 +0.16 0.03 +0.11 0.07 +0.03 0.03 +0.14 0.06 +0.05 0.02 +0.05 0.03 2.9 ITER 

HIP 55486 5371 44 4.51 0.06 +0.46 0.03 +0.42 0.03 +0.58 0.07 +0.37 0.05 +0.39 0.02 +0.48 0.03 9.6 ITER 

HIP 56729 5421 44 4.47 0.06 -0.09 0.03 -0.04 0.03 -0.10 0.03 -0.02 0.05 -0.04 0.02 -0.11 0.03 2.7 ITER 

HIP 56930 5025 44 4.61 0.06 -0.12 0.03 -0.17 0.05 -0.20 0.03 -0.15 0.05 -0.14 0.02 -0.18 0.03 7.0 ITER 

HIP 59080 5423 44 4.41 0.06 -0.16 0.03 -0.13 0.03 -0.12 0.03 -0.16 0.05 -0.11 0.02 -0.16 0.03 4.9 ITER 

HIP 61081 5298 44 4.59 0.06 -0.54 0.03 -0.42 0.03 -0.49 0.03 -0.33 0.05 -0.26 0.02 -0.48 0.03 2.5 ITER 

HIP 61189 4655 44 4.93 0.06 +0.11 0.03 +0.05 0.03 -0.01 0.03 +0.02 0.05 +0.08 0.02 +0.09 0.03 16.6 VESTA 

HIP 61589 5700 44 4.49 0.06 -0.05 0.03 -0.08 0.05 -0.08 0.03 -0.05 0.05 -0.05 0.02 -0.06 0.03 2.3 ITER 

HIP 64345 5495 44 4.40 0.06 -0.57 0.03 -0.39 0.03 -0.49 0.03 -0.33 0.05 -0.25 0.02 -0.52 0.03 1.8 ITER 

HIP 64797 5041 44 4.60 0.06 -0.12 0.03 -0.13 0.03 -0.16 0.03 -0.11 0.05 -0.16 0.02 -0.17 0.03 15.1 ITER 

HIP 65636 4619 44 4.64 0.06 +0.13 0.03 +0.01 0.03 -0.11 0.03 -0.08 0.05 +0.00 0.02 +0.02 0.03 13.4 ITER 

HIP 65963 5478 44 4.52 0.06 -0.14 0.03 -0.08 0.03 -0.19 0.03 -0.08 0.05 -0.09 0.02 -0.14 0.03 1.7 ITER 

HIP 68799 5492 44 4.42 0.06 -0.03 0.03 -0.06 0.03 +0.00 0.03 -0.05 0.05 -0.02 0.02 -0.04 0.03 2.7 ITER 

HIP 70100 4886 44 4.57 0.06 +0.15 0.03 +0.13 0.03 +0.26 0.03 +0.09 0.05 +0.12 0.02 +0.16 0.03 9.9 ITER 

HIP 70426 4801 44 4.58 0.06 +0.09 0.03 +0.03 0.04 +0.13 0.04 -0.00 0.05 +0.06 0.02 +0.10 0.03 12.8 ITER 

HIP 74396 5124 44 4.60 0.06 -0.09 0.03 -0.07 0.03 -0.08 0.03 -0.06 0.05 -0.12 0.02 -0.13 0.03 3.5 ITER 

HIP 74734 5822 44 4.34 0.06 -0.32 0.03 -0.27 0.03 -0.37 0.03 -0.20 0.05 -0.23 0.02 -0.39 0.03 2.1 ITER 

HIP 75069 5196 44 4.61 0.06 -0.38 0.03 -0.34 0.03 -0.36 0.03 -0.33 0.05 -0.33 0.02 -0.40 0.03 2.8 ITER 

HIP 76668 4636 44 4.65 0.06 -0.06 0.03 -0.14 0.03 -0.22 0.06 -0.17 0.05 -0.10 0.02 -0.12 0.03 13.7 ITER 

HIP 78969 4899 44 4.58 0.06 +0.19 0.03 +0.12 0.03 +0.23 0.06 +0.08 0.05 +0.13 0.02 +0.20 0.03 9.0 ITER 

HIP 79629 5600 50 4.48 0.06 -0.25 0.04 -0.25 0.06 -0.28 0.03 -0.25 0.05 -0.20 0.02 -0.31 0.03 2.0 ITER 

HIP 84616 4826 44 4.62 0.06 -0.12 0.03 -0.09 0.03 -0.09 0.03 -0.09 0.05 -0.10 0.02 -0.14 0.03 6.4 ITER 

PM J1742+1645 5492 44 4.50 0.06 -0.09 0.03 +0.01 0.03 -0.06 0.03 -0.00 0.05 -0.11 0.02 -0.15 0.03 4.2 VESTA 

HIP 87082 5795 44 4.40 0.06 -0.05 0.03 -0.03 0.05 -0.17 0.03 +0.02 0.05 -0.02 0.02 -0.09 0.03 2.7 ITER 

HIP 88188 5299 53 4.54 0.06 +0.05 0.04 -0.04 0.07 -0.06 0.03 -0.05 0.13 +0.02 0.03 +0.04 0.03 3.2 ITER 

HIP 88365 5371 44 4.59 0.06 -0.66 0.03 -0.34 0.03 -0.52 0.03 -0.33 0.05 -0.33 0.02 -0.57 0.03 3.7 ITER 

HIP 90246 4494 44 4.69 0.06 -0.01 0.03 -0.08 0.03 -0.26 0.05 -0.08 0.05 -0.08 0.05 -0.09 0.03 11.3 ITER 

HIP 104097 4428 44 4.72 0.06 -0.38 0.03 -0.27 0.03 -0.51 0.03 -0.15 0.05 -0.37 0.04 -0.40 0.03 15.1 ITER 

PM J2206+4322E 5652 44 4.62 0.06 +0.30 0.03 +0.15 0.03 +0.37 0.03 +0.22 0.06 +0.29 0.02 +0.33 0.03 6.2 VESTA 

HIP 114156 4314 66 4.72 0.06 -0.04 0.08 -0.06 0.05 -0.37 0.15 -0.08 0.05 -0.24 0.07 -0.15 0.03 15.5 ITER 

HIP 117960 5252 44 4.54 0.06 +0.12 0.03 +0.06 0.03 +0.14 0.03 +0.04 0.05 +0.12 0.02 +0.11 0.03 4.0 ITER 

HIP 118282 5178 44 4.64 0.06 -0.67 0.03 -0.47 0.03 -0.54 0.03 -0.36 0.05 -0.35 0.03 -0.57 0.03 2.0 ITER 

a ITER: parameters determined using Hipparcos parallaxes and Y 2 isochrones. VE STA: parameters determ ined using classical SME 
fitting (no parallax information included) with a correction using Vesta as described in |Valenti &: F ischer (2005]|. 

of the control sample, (3) the mean difference between 
the SPOCS metallicities and the literature metallicities 
in the control sample A[Fe/H] control < 0.07 dex, and 
(4) the resulting scatter in the control sample cr contro i 
< 0.08 dex. These limits are designed to keep uncertain- 
ties in the primary star metallicities well below the pre- 
cision already obtained for determining M dwarf metal- 
licities (e.g., R12 and T12). We adopt <7 con troi as the 
uncertainty in [Fe/H] for a given literature source. As 
an example , we show metallicities from both f~ 



of metallicities, the adopted systematic offset (which we 
apply for all calculations in this paper) for that source, 
and the adopted uncertainty in Table [3] 

Late-K/M Dwarf Spectral Types 
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et al.| ( |2007| ) and SPOCS in Figure [I] We list all sources tral types accurate to 



We determine M dwarf spectral types using indices 
at both visible and NIR wavelengths. We use t he em- 
pirica l spec tral type-band index relations from |Lepine| 
et al. (2012), whi ch have been calibra ted to work on the 
]STOFS/TJH2.2m. 



Lepine et al 
UTTu 



(2012) determined spec- 



Dtypes based on empir- 
ical relations between spectral type and the strengths 
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Figure 1. Comparison of metallicities from Ramirez et al.] 1 2007[ l 
and those from SPOCS ( |VaIenti k. Fischer|20u5] | tor the 112 stars 
present in both samples! Ine [f 'e/rl] values are quite consistent 
once corrected for a systematic offset of 0.08 dex. The remaining 
scatter is only ~ 0.06 dex. Assuming SPOCS [Fe/H] v alues are 
accurate to .03 dex, this implies [Fe/H] metallicities from |Ranu'rez| 
|et al.| ( |2007| are accurate to ~ 0.05 dex. We perform a similar 
analysis lor all other metallicity sources (see Table 



Table 3 

Corrections Applied to Primary Star Metallicities 



Source 


Control Stars 


A [Fe/H] 


•^control 


No. of Stars Used" 


VF05 




0.00 


0.03 


33 


CFHT 




0.00 


0.03 6 


50 


C01 


294 


0.02 


0.07 


1 


M04 


178 


0.04 


0.08 


1 


LH05 


174 


0.02 


0.06 


1 


T05 


127 


0.01 


0.05 


1 


B06 


33 


0.07 


0.07 


1 


Ra07 


112 


0.08 


0.05 


1 


Ro07 


127 


0.00 


0.07 


5 


F08 


165 


0.03 


0.06 


2 


Cll c 


614 


0.00 


0.08 


10 


Sll 


50 


0.03 


0.04 


1 


N12 d 


125 


0.00 


0.05 


2 



See Table [T] for the abbreviations on references. "Number of stars 
from listed source used in our final wide binary sample. 
'Typical uncertainty. Errors for individual stars listed in Table [2] 
c Metallicities determined from S tromgren photometr y. 
d F3ased on a control sample from |Santos et ah] j2004"| l and |Sousa | 
|et al.||20TT| |, on which abundances in IN 12 are anchored. 



of TiO and CaH bands (Reid et al. 1995). It has been 
shown that CaH is sensitive to spectral type, and that 
TiO is sensitive to both sp ectral type and metallicity 
(Woolf & Wallerstein 2006). As a result, for stars with 
[Fe/H|< —0.5 (the metallicity of the primary star) we 
base our visible wavelength spectral types solely on rela- 
tions using CaH bands. 

R12 showed that one can determine temperatures and 
spectral t ypes using a mod ified version of the H 2 0-K 
(H 2 Q-K2, [Covey et al.|2 010) index. R12 calibrated their 
spectral types based on A -band spectra of stars from the 
Resear ch Consortium on Nearby Stars Measuring (RE- 
CONS; [Henry et al.||1994[ ). Their calibration is accurate 
to ~ 0.6 subtypes. 

Figure [2] compares the spectral types derived from vis- 
ible wavelength indices versus those derived using the 
H2O-K2 index. Although there is good agreement be- 
tween the two techniques for the later-type stars in our 
sample, for stars earlier than Ml (as determined by TiO 



H 



GO 



M6 

M4 
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M0 
K5 



+ 



K5 M0 M2 M4 M6 

Visible Wavelength Spectral Type 

Figure 2. Spectral types derived from the H2O-K2 index (RA12) 
vs. those derived from TiO and CaH indices at visible wavelengths. 
The dashed line indicates a perfect agreement and typical errors 
are shown in the bottom right of the plot. For spectral types later 
than Ml the different techniques agree, but for earlier-type stars 
the H2O-K2 index tends to return a later spectral type than visible 
wavelength indices. 

and CaH bands), spectral types determined from H2O- 
K2 are systematically later than those from visible wave- 
length indices. The H2O features become quite weak in 
the spectra of late-K and early-M stars and R12 caution 
using it on stars with T e g > 4000. Further, the spectral- 
type calibration from R12 does not include any K stars, 
and is therefore unreliable for the warmest stars in our 
sample. As a result, we choose to use spectral types 
determined from our visible wavelength spectra for the 
entire sample. 

5. IDENTIFYING METAL-SENSITIVE INDICES 

To determine which features in the companion dwarf 
spectra best correlate with metallicity we perform a sys- 
tematic analysis of our sample of spectra and metallici- 
ties. Our analysis proceeds as follows: 

1. A center wavelength is selected, starting at the 
blue end of the spectrum (~ 0.33 /mi) and incre- 
mentally increasing by 0.00015 /mi (1.5 A) after all 
other steps are complete. This process is repeated 
until the center is at the red end of the spectrum 
(~ 2.4 /im) and excludes the gap in all SpeX spec- 
tra at 1.85 /im. 

2. For each feature center, we select a feature width 
starting at 0.002 /mi (20 A), and then increased in- 
crementally by 0.00015 /im (1.5 A) after completing 
all following steps. We use 20 A as a minimum, as 
features smaller than this have considerable Pois- 
son noise (making their measurement difficult). We 
use an upper limit of 0.01 /im (100 A) for the fea- 
ture width, as regions of the spectra larger than 
this likely contain multiple features that should be 
treated separately. 

3. The equivalent width is calculated for each feature 
using the approximation: 



EW> 



E 



1 



F c (Xi) 



AAi, 



(3) 



where Xi is the wavelength at pixel i, F is the flux 
at Xi, F c is the pseudo-continuum at Xi, and the 



Prospecting in late-type dwarfs 



7 



sum is computed over all n pixels within a given 
feature. We compensate for the low resolution by 
interpolating the spectrum near the edge to a much 
higher resolution (R > 10,000). We experiment 
with different techniques to calculate the pseudo- 
continuum (see below). The list of continuum re- 
gions used is listed in Table |4j many of which are 
taken from T12. 

4. A temperature-sensitive parameter r, is calculated 
from each spectrum, r is defined based on the cen- 
ter wavelength of the selected feature from step (2). 
Specifically, r is set to be the if-band H2O (H 2 0- 
K2) index as defined by R12 if the feature is cen- 
tered in the K band, the H band H 2 (H 2 0-H) 
index as defined by T12 if th e feature is in the H 
band, or the Color 1 index from Hawley et al.| ( 2002) 
if the feature is drawn from visible wavelengths. If 
the feature falls within the J-band, we use a new 
H2O-J index defined as: 



Table 4 

Continuum Regions Used 



H 2 0-J 



(JT(1.210- 1.230))/(J"(1.313- 1.333)) 
(J"(1.313 - 1.313))/(J"(1.311 - 1.351)) ' 



(4) 



where (^(a—b)) indicates the median flux level in a 
wavelength range between a and b (in /mi). H 2 0-J 
is defined to select regions relatively clear of atomic 
or molecular features and to correlate well with the 
H 2 0-K and H 2 0-H indices. 

5. Using least-squares, the best fit is found for the 
equation: 



[Fc/H]. i = A + B x EW, + Cx n 



(5) 



where [Fe/H]i is the metallicity of ith primary star 
(assumed to be the metallicity of the companion M 
dwarf), EWi is the calculated equivalent width of 
the selected feature in the ith late-K or M dwarf 
companion spectrum, and A, B, and C are fitting 
parameters. The quality of the fit is measured by 
the adjusted square of the multiple correlation co- 
efficient (i?a P )j which is defined as: 



7->2 _ -I _ i n ~ !)E(^ 

(n- P )J2(y 



1 . model Vi ) 



(6) 



where p is the number of changeable parameters 
(i.e. A, B, and C), n is the number of data points in 
the fit, yi is the metallicity of the ith primary star, 
Vi, model is the metallicity of the ith star predicted 
by the fit, and y is the average of yi. A R 2 p closer 
to 1 implies that the model accurately explains 
the variance of the sample, while R ap =0 implies 
that it can explain none. For Equation p = 3, 
yj=[Fe/H]j, and yi, mo dei = A + B x EW, +Cxt,. 
Note that for n > p, R 2 ~ R 2 . 



To asses the significance of the assigned R ap value, 
[Fe/H] (or [M/H]) values are randomly reassigned 
among the stars, and step 5 is repeated 1000 times 
(re-randomizing the metallicities each time). The 
resulting distribution of the 1000 R ap values gives 
the level above which the R 2 value (determined 



Visible 


J-Band 


_ff-Band 


iC-Band 


W 


(Mm) 


(H 


(cm) 



0.4035- 
0.4135- 
0.4425- 
0.4610- 
0.4680- 
0.5269- 
0.5660- 
0.6586- 
0.7041- 
0.7390- 
0.8100- 
0.8230- 
0.8590- 
0.8890- 
0.9100- 
0.9220- 



0.4080 
0.4180 
0.4450 
0.4625 
0.4700 
0.5299 
0.5675 
0.6607 
0.7049 
0.7500 
0.8160 
0.8300 
0.8620 
0.8920 
0.9120 
0.9255 



0.9790- 
1.0610- 
1.1260- 
1.1530- 
1.1890- 
1.2140- 
1.2250- 
1.2550- 
1.2700- 
1.2950- 
1.3040- 
1.3214- 
1.4090- 



0.9890 
1.0650 
1.1300 
1.1580 
1.1930 
1.2180 
1.2300 
1.2634 
1.2730 
1.2970 
1.3070 
1.3270 
1.4150 



1.4440- 
1.4644- 
1.4921- 
1.5060- 
1.5190- 
1.5920- 
1.6230- 
1.6935- 
1.7530- 



1.4480 
1.4710 
1.4965 
1.5090 
1.5220 
1.5960 
1.6310 
1.6980 
1.7570 



1.8860- 
1.9350- 
1.9610- 
2.0500- 
2.0800- 
2.1330- 
2.1530- 
2.1670- 
2.1940- 
2.2130- 
2.2450- 
2.2717- 
2.2850- 
2.3050- 
2.3600- 
2.3710- 
2.3950- 



1.8900 
1.9400 
1.9700 
2.0540 
2.0870 
2.1351 
2.1590 
2.1720 
2.1985 
2.2190 
2.2520 
2.2781 
2.2900 
2.3105 
2.3640 
2.3760 
2.4050 



from non-random metallicities) can be considered 
significant. We consider the given feature center 
to be a bona-fide metal sensitive feature if R^ p is 
higher than the 99.9% highest R^ p value from the 
randomly assigned metallicities (henceforth R^ and ). 



For each increment in feature center and width, we record 
the resulting R ap and R^nd- We show the resulting dis- 
tribution of R 2 and R 2 and values as a function of the 
feature's central wavelength in Figure [3j 

We repeat our analysis using various methods of fit- 
ting for the pseudo-continuum. Better estimates of the 
continuum should result in more accurate line measure- 
ments, and therefore higher R ap values for the same fea- 
tures. In one experiment we tried to fit the global spec- 
trum with a high order (> 10) polynomial. We tested 
fitting each band (visible, JHK) with 3rd through 6th 
order polynomials, as well as with 3rd through 6th order 
Legendre polynomials. Interestingly, we found we had 
the best (highest R ap with respect to R 2 and ) results when 
fitting the pseudo-continuum using a linear fit of the con- 
tinuum regions immediately blueward and redward of the 
selected feature. We use this fitting procedure for all cal- 
ibrations derived in Section |6] 

We run additional experiments to test the influence 
of spectral type on determination of metallicities: we 
repeat our analysis on just the early-type stars in our 
sample (K5.5-M2.0) and again with just the late-type 
stars (M2.0-M6). The split roughly corresponds to our 
median spectral type (~ M2). It is possible to parse 
our data into smaller spectral type ranges, although this 
will proportionately shrink each sample. Thus this would 
make it difficult to identify metal sensitive features that 
have not been previously discovered. 

Lastly, we rerun our analysis using [M/H] rather than 
[Fc/H] for the metallicity of the primary star. R12 were 
able to derive better fits between if-band atomic lines 
and [M/H] than for [Fe/H], perhaps due to variations 
in [a/Fe] creating discrepancies between the actual mea- 
sured quantities (Na and Ca in the case of R12) and Fe 
abundances. However, many of the literature sources we 
draw from do not provide [M/H] values for the primary 
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Figure 3. Results of our systematic search to find the features that best correlate with [Fe/H] in late-K or M dwarf spectra. Feature 
centers (shown on the X axis) and widths are changed incrementally covering a range of widths (20 A to 100 A) and feature centers (0.35 /im 
to 2.4 /mi). Red points indicate the Rj? and values derived from randomly reassigning primary star metallicities and repeating our process 
1000 times (the 99.9% highest resulting R ap values are shown). In this particular analysis, we fit for [Fe/H], and fit the pseudo-continuum 
using just the continuum regions immediately blueward and redward of a given feature. A range of feature widths are shown, which results 
in some range to the distribution for a given central wavelength. Although there are > 400,000 unique combinations of feature center and 
width, only 15,000 are shown for simplicity. 



stars (just [Fe/H]). As a result, there are 5 fewer binaries 
when using [M/H] and 16 primary stars have a different 
source for their [M/H] values than their [Fe/H] values 
(see Table [l] for a full list of [Fe/H] and [M/H] sources 
used) . 

We consider a feature to be metal-sensitive if a given 
feature's center and width has an R^ p value are above its 
corresponding Rj 2 and value. This criterion ensures that 
no features are identified simply by coincidence. 

Once features are identified, we then attempt to derive 
a calibration for a given wavelength range (e.g., J-band) 
by solving the equation (by least squares): 



[Fe/H] j = A 



N 



(7) 



where [Fe/H] j is the metallicity of the ith primary star, N 
is the total number of features of interest among M\ fea- 
tures identified as metal-sensitive in a given wavelength 
range, Tj is the temperature sensitive parameter selected 
based on the wavelength regime (see above), EWjj is the 
equivalent width of the jth feature measured for the ith 
star, and A, C, and the Bj's are fitting parameters. We 
find the best fit for iV=l,2,3,...< M\ until the increase 



in R ap is negligible (Ai? ap < 0.03) or it is clear from 
visual inspection of the data that the adding of further 
variables is over fitting the data. This limit is usually hit 

at TV = 3 4. Although we use [Fe/H] in Equations [5] 

andJTJ we also perform the same procedure using [M/HJ7 

6. DETERMINATION OF M DWARF METALLICITIES 

The first thing our analysis gives us is a catalog of 
metal-sensitive features, which we list in order of A c in 
Table [5] In total we find 120 features that are statisti- 
cally significant predictors of metallicity, although only 
20 of these are used in our final calibrations. We iden- 
tify a number of previously known metal-sensitive fea- 
tures, as well as many of new ones. One of the most 
metal-sensitive features is the Na I doublet in the K- 
band (2.208 /an), already identified by R10. Our analy- 
sis identifies the Ca I (1.616 /imand 1.621 /mi) and K I 
(1.5176 /im) lines shown to be metal-sensitive by T12. 
In fact, our analysis locks on to very similar wavelength 
centers and widths as those found by T12 for both H- 
band and if-band features. Since our analysis covers 
all wavelengths and is completely blind (e.g., they have 
no a-priori line lists or knowledge of feature size) this 
suggests that our purely empirical analysis is identifying 
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Table 5 

Metal-sensitive Features 





Center A 


Width 






R 


2 

IP 








[fan] 


[A] 




[Fe/H] 




[M/H] 










All 


Early 


Late 


All 


Early 


Late 


F01 


0.4648 


23 




0.56 






0.54 




F02 


0.5608 


20 




0.43 






0.30 




F03 


0.6118 


20 


0.14 


0.20 


0.26 


0.15 


0.26 


0.39 


F04 


0.6232 


20 






0.28 






0.41 


F05 


0.6416 


41 


0.13 




0.42 


0.22 




0.57 


F06 


0.7540 


20 


0.30 


0.51 


0.39 


0.28 


0.44 


0.38 


F07 


0.8208 


35 


0.36 


0.77 


0.23 


0.31 


0.69 


0.22 


F08 


0.8684 


26 


0.14 


0.63 


0.23 


0.15 


0.53 


0.18 


F09 


1.1396 


26 


0.35 


0.46 


0.28 


0.27 


0.35 


0.19 


F10 


1.2698 


98 


0.42 


0.47 


0.37 


0.43 


0.47 


0.35 


Fll 


1.2908 


20 


0.51 


0.63 


0.41 


0.41 


0.55 


0.26 


F12 


1.3148 


50 


0.39 


0.72 


0.23 


0.31 


0.66 


0.25 


F13 


1.3344 


23 


0.15 


0.42 


0.31 


0.16 


0.35 


0.38 


F14 


1.4766 


41 


0.54 


0.60 


0.44 


0.45 


0.51 


0.29 


F15 


1.4836 


23 


0.43 


0.73 


0.18 


0.36 


0.66 


0.25 


F16 


1.5172 


33 


0.48 


0.70 


0.48 


0.48 


0.69 


0.45 


F17 


1.6158 


23 


0.60 


0.88 


0.28 


0.55 


0.86 


0.22 


F18 


1.7261 


32 


0.24 


0.39 


0.26 


0.18 


0.38 


0.19 


F19 


2.2079 


68 


0.86 


0.88 


0.83 


0.78 


0.72 


0.68 


F20 


2.3242 


38 


0.63 


0.75 


0.58 


0.56 


0.57 


0.50 


F21 


2.3342 


35 


0.61 


0.73 


0.55 


0.57 


0.58 


0.49 


F22 


2.3844 


35 


0.54 


0.69 


0.35 


0.50 


0.51 


0.32 



Here we show only the features used in our final calibrations. The 
full version of the table with all metal-sensitive features identified 
by our analysis will be available electronically. 

a ... denotes that the feature did not have an R ap value above the 
R^ an , value, and thus is not considered a statistically significant 
metal-sensitive feature. 

b Full: K5.5-M6.0, early: K5.5-M2.0, late: M2.0-M6.0. 



metal-sensitive atomic and molecular lines. 

Because our only restriction was the size of the fea- 
ture (< 100 A), our technique can easily identify areas 
of the spectrum corresponding to several, even unrelated 
lines. Some of these regions may be associated with dou- 
blets/triplets from the same atomic species, with broad 
molecular bands, or with sets of lines that are blended 
at our resolution. Some features in Table [5] may not 
correspond to any one specific element or molecule, but 
simply to a region of the spectrum that undergoes overall 
changes as a function of the metallicity of the star. 

We show the distribution of R 2 p values as a function 
of feature width and center for two example wavelength 
regions in Figure [4] Interestingly, features in the H and 
K band features yield similar R 2 p for a wide range of fea- 
ture widths. However, the opposite is seen in the visible 
end of the spectrum, where features perform better near 
the minimum feature width (20 A). This is most likely 
due to crowding at visible wavelengths. 

We find far better metal- sensitive features by doing a 
separate analysis for earlier- type dwarfs (K5.5-M2) and 
for later type dwarfs (M2-M6). Specifically, the best fit 
we achieve when fitting for [Fe/H] using all stars in our 
sample yields R^p=0.54, whereas when we split up the 
sample by spectral type we achieve R 2 p =0.84 for K5.5- 
M2.0 and R 2 p =0.68 for M2.0-M6.0. This is not unex- 
pected, many of the most metal-sensitive features for 
K5.5-M2 dwarfs become blended with molecular bands 
(which grow as a function of spectral type) at the res- 
olution of SNIFS. Further, features blueward of 0.5 fim 
tend to have very low S/N for stars later than M2, and 
are not to be very useful. This also suggests that bet- 
ter results could be achieved on later type M dwarfs with 



modest improvements in resolution, to better distinguish 
the lines. 

We find the best empirical fits to Equation [7] for each 
wavelength regime (visible, J, H, and K) and metalli- 
cally metric ([Fe/H] and [M/H]). They are: 



[Fe/H]y, e 


= 0.53^7 + 0.26^1 - 0.16F 02 


(8) 




-0.784(Colorl) - 0.34 




[M/H] v ,c 


= 0.38^07 + 0.21^1 + 0.29F 08 


(9) 




-0.504(Colorl) - 0.79 




[Fe/H] v ,i 


= -0.20F u5 + 0.48F 08 + 0.24F 07 


(10) 




+0.14^3 - 0.204(Colorl) - 0.32 




[M/H] v ,i 


= -0.065F 05 - 0.071F 04 - 0.30Fo6 


(11) 




+0.719(Colorl) - 0.24 




[Fe/H] j 


= 0.29Fi + 0.2LF g + O.26F12 


(12) 




-O.26F13 - 0.190(H 2 O-J) - 1.03 




[M/H] j 


= 0.32F 10 + 0.46^11 + 0.076F 09 


(13) 




+1.213(H 2 0-J) - 1.97 




[Fe/H]fj 


= 0.40Fi 7 + 0.51^14 - 0.28F 18 


(14) 




-1.460(H 2 O-H) + 0.71 




[M/H] H 


= 0.38Fi 7 + 0.40^16 + O.4IF15 


(15) 




+0.194(H 2 O-H) -0.76 




[Fe/H] K 


= 0.19Fi 9 + 0.069F 22 + 0.083F 20 


(16) 




+0.218(H 2 O-K) - 1.55 




[M/H] K 


= 0.12Fi 9 + 0.086F 22 + 0.13F 2 i 


(17) 




+0.245(H 2 O-K) - 1.18 





where F=#= refer to features listed in Table [5j the sub- 
scripts refer to the wavelength bands where the calibra- 
tion is useful (V referring to visible wavelengths). An 
additional subscript is added (e or 1) for calibrations in 
visible wavelengths to denote which formula is valid for 
early (K5.5 to M2) and late (M2 to M6) dwarfs. We show 
the primary star metallicity as a function of the derived 
metallicity for the companion dwarf for each of these 10 
calibrations in Figure pi and list reduced \ 2 j R-ap > r0 °t 
mean square error (RMSE), in Table [6] The RMSE indi- 
cates how useful a model is at prediction (lower numbers 
indicate the fit is a better predictor) and is defined as: 



RMSE = 



\ 



(Vi, 



model 



(n - p) 



(18) 



Lower R 2 p and higher RMSE values may in part be 
due to differences in S /N as a function of wavelength. We 
estimate measurement noise sources by adding synthetic 
noise to each spectrum consistent with the observed S /N 
of that spectrum, then recalculating the metallicity of 
the M dwarf using the appropriate equation above. The 
standard deviation in the metallicity estimate from 1000 
different additions of noise pattern is assumed to be the 
measurement error. This error is what is what we use 
for our calculation of the reduced x 2 f° r each fit. Thus 
the reduced \ 2 values probe how much of the noise comes 
from measurement (errors)A reduced x 2 close to 1 would 
suggest that most or all of the error from the fit is due 
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Figure 4. Distribution of R^ p values as a function of feature center and width for the iC-band (left) and red end of the visible wavelengths 
(right). The visible wavelength distribution is based on our analysis using just K5.5-M2 dwarfs. Data points which do not have R^ p >R^ and 
are white (i.e. not shown). Surprisingly, there is not much change in R ap as a function of the selected width for a given feature in the 
A"-band, however we do see a trend towards smaller widths in the visible wavelengths. 



to measurement (e.g. Poisson) noise. 

The dependence on r varies significantly between equa- 
tions. This is most likely due to different features cap- 
turing some of the temperature-dependence and/or that 
the Colorl and H2O indices are not accurately model- 
ing temperature dependencies across the full sample. It 
is also interesting that some coefficients of features are 
negative. This may be due to a combination of factors, 
including changes in [Fe/H] versus [a/Fe] (e.g. Ca is an 
a element) or complex relations between T e g and [Fe/H], 
that vary for each feature. Whatever the case, since these 
fits are purely empirical, we should be cautious not to 
over-interpret the physical meaning of any particular co- 
efficient or feature. 



7. 



ASSESSING AND RECALIBRATING EXISTING 
TECHNIQUES 



In addition to defining our own metallicity calibrations, 
we can use our sample to test existing metallicity esti- 
mators, as well as improve the existing calibrations. Like 
before, we use R ap , andRMSE as our standard metrics 
to asses the quality of a calibration. We summarize our 
refits in Table 

7.1. CTiO/CaH 

Much effort has gone into determine M dwarf metallic- 
ities using visible wavelength spectra. Most of this has 





Lepine et al. 


2007 


Woolf et al. 


2009| Dhital et al. 


2012 


. However, tl 


le setup of our ana 


ysis means that we 



on spectroscopic indices (not equivalent widths). Band 
indices (e.g., Ti05, CaH3, etc.) are calculated from the 
ratio of the flux in region a to the flux in region b using 
the approximation: 



Index 



[Ei= fe -F(Ai)]/W 



(19) 



where w a and Wb are the widths of region a and b in 
angstroms, respectively. The sums are computed over 
all pixels i in region a and b, respectively. Our analysis 
only makes use of equivalent widths (see Equation [3]) . 
Further, we do not allow high order terms, while £ gen- 
erally requires 3rd or 4th order polynomials of the CaH 



index (e.g. Lepine et al. 2012). However, this does not 
prevent us from using our data to test the performance 
of C- 

We calculate the CaH2, CaH3, and Ti05 indices fol- 



lowing the definitions from Reid et al. (1995). We com- 
pute corrected in dices (CaH2 e CaH3 e , and li05 e ) using 
the formula from Lepine et al. (2012), which include cor- 
rections for the SNIFS instrument. We use these to com- 
pute C following the formula as defined by |Lepine et al.| 
p07| : 

1 - Ti05 

C = l-[Ti05V (20) 

where [Ti05] z@ is a function of CaH = CaH2 + CaH3. 
We use the formula for [Ti05]z o 
p0l2|. 



from Lepine et al. 



[TiO5] Z0 =0.622 - 1.906(CaH c ) 

+2.211(CaH c ) 2 - 0.588(CaH c ) 3 . 



(21) 



We plot the primary star metallicities as a function of 
the derived ( values as filled points in Figure [6] ( shows 
a weak trend with metallicity in both [Fe/H] and [M/H]. 
From this we derive the following relationships: 



[Fe/H] : 
[M/H] : 



:0.98C 
:0.68C 



1.04 
0.74 



(22) 
(23) 



Like before, we randomly reassign the metallicities to 
different CPM pairs, and attempt to compute an R^ an( j 
value. We find that both Equations [23] and [23] give 
R ap <R^ and . Further, the ( parameter only correctly 
identifies one companion as a subdwarf (LHS 1812/PM 
I06032+1921S). Although this is the most metal-poor 
star in our sample, there are 12 other stars in our sam- 
ple with [Fe/H]< —0.5 but ( values c onsistent with sola r 
metallicity « > 0.825). According to |Woolf et ak] ( |2009[ ), 
M dwarfs with [Fe/H]< -0.34 should have C <TI825, 
and be labeled as sdM, suggesting a problem with these 
stars. If we remove these pairs, we derive the following 
relations: 



[Fe/H] = 1.26C - 1.25 
[M/H] =0.88C- 0.89 



(24) 
(25) 
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M Dwarf Metallicity [Fe/H] M Dwarf Metallicity [M/H] 



Figure 5. Metallicity for the primary star vs. derived metallicity of the late K or M dwarf companion based on Equations (T8P— ( fXTP - 
Calibrations for [Fe/Hl are shown on the left plots, while those for [M/H] are shown on the right plots. Statistics on the quali ty oT the tit 
can be found in Table Y axis error bars shown are based on la Gaussian errors for the primary star metallicity (see Section |4.1| . Error 
bars for the K/M dwarf metallicity are the la standard deviation of 1000 recalculations of the K/M dwarf metallicity after adding noise to 
each spectrum consistent with its S/N. 



Table 6 

Metallicity Calibration Statistics 



Equation # Band 


Wavelength Range 


SpT Range 


Metal Type 


R-ap 


RMSE 


a 




^Rod 


8 




Optical 


0.35< A < 1.00 


K5.0-M2.0 


[Fe/H] 


0.84 


0.07 





13 


8.8 


9 




Optical 


0.35< A < 1.00 


K5.0-M2.0 


[M/H] 


0.80 


0.06 





11 


6.4 


T 


1 


Optical 


0.35< A < 1.00 


M2.0-M6 


[Fe/H] 


0.68 


0.06 





14 


7.1 


11 


Optical 


0.35< A < 1.00 


M2.0-M6 


[M/H] 


0.65 


0.06 





11 


7.8 


12 


J 


1.00< A < 1.44 


K5.0-M5.0 


[Fe/H] 


0.71 


0.07 





16 


11.4 


H 


J 


1.00< A < 1.44 


K5.0-M5.0 


[M/H] 


0.55 


0.08 





15 


9.4 


11 


H 


1.44< A < 1.80 


K5.0-M5.0 


[Fe/H] 


0.77 


0.07 





14 


3.7 


To 


H 


1.44< A < 1.80 


K5.0-M5.0 


[M/H] 


0.73 


0.06 





12 


4.2 


16 


K 


1.80< A < 2.45 


K5.0-M5.0 


[Fe/H] 


0.86 


0.06 





11 


1.7 


17 


K 


1.80< A < 2.45 


K5.0-M5.0 


[M/H] 


0.77 


0.05 





10 


3.8 



Table 7 

Assessment of Previous Metallicity Indicators 



Technique 
CTiO/CaH 
CTiO/CaH 

J — K 

J -K 

//-Band 

//-Band 

K-Band 

ii-Band 



Type 

[Fe/H] 

[M/H] 

[Fe/H] 

[M/H] 

[Fe/H] 

[M/H] 

[Fe/H] 

[M/H] 



0.58 
0.61 
0.30 
0.25 
0.74 
0.71 
0.76 
0.75 



RMSE 

0.28 

0.23 

0.19 

0.16 

0.14 

0.12 

0.14 

0.13 



Eqn No. 

W 

27 

29 
30 
31 
32 
33 
31 



These formulae are highly significant; they yield R^ p val- 
ues of 0.58 and 0.52, and RMSE values of 0.22 and 0.20, 
respectively. This suggests that ( may be useful at pre- 
dicting metallicitie s for [Fe/H]> +0.05 (the limit of the 



Woolf et al. (2009) calibration), provided the high-zeta 
lo w-metaUicity stars c an be explained. 

Woolf et al] ( |2009| using R ~ 3000 spectra) derive 
a relation between metallicity and the £ parameter us- 
ing a mix of wide binaries (for which the primary star 
metallicity is known) and single (K/M) stars with high- 
resolution spe ctra, analyzed using the MOOG software 
dSnedenl 19731) with NEXTGEN models (IHauschildt et al. 
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TableJ 

Stars Observed from |Woolf et al.| ( [2009] ) 



IN eLIIlL. 


[rc/nj 




a 

[Fe/H] 


[ivi/nj 




LHS 38 


-0.43 





05 


-0.40 


0.94 


HIP 1386 


0.16 





10 


0.15 


0.98 


HIP 59514 


-0.05 





08 


-0.03 


1.06 


HIP 89490 


-0.53 





08 


-0.44 


0.85 


HIP 98906 


-0.62 





10 


-0.52 


0.56 


HIP 105932 


-0.37 





05 


-0.30 


0.87 


HD 18143B 


0.19 





11 


0.18 


1.03 


HD 88230 


-0.03 





18 


-0.05 


1.00 


HD 95735 


-0.42 





07 


-0.40 


1.10 


GJ 129 


-1.66 





05 


-1.33 


0.08 


GJ 1177B 


-0.09 





10 


-0.06 


0.91 


GJ 3212 


-0.08 





05 


-0.06 


0.70 


-Lxlo O04 


1 A 1 

-1.41 





A 1 
01 


-1.15 


0.00 


GJ 9722 


-0.83 





04 


-0.70 


0.58 


LHS 174 


-1.11 





05 


-0.95 


0.57 


LHS 182 


-2.15 





03 


-1.88 


-0.27 


LHS 491 


-0.93 





08 


-0.78 


0.45 


LHS 3084 


-0.73 





05 


-0.64 


0.78 


LHS 156 


-1.00 





04 


-0.85 


0.60 


LHS 161 


-1.30 





04 


-1.06 


0.29 


LHS 318 


-1.26 





05 


-1.03 


0.25 


LSPMJ2205+5353 


-1.29 





08 


-1.06 


0.12 



|stein| ( |200e] ) , or[Wooll ct al. | |2U09[ |. 

"Determined from SNIFS spectra. 



Waller-| 



19991. Although many of their wide binaries are also in 
our binary sample (including LHS 1812), there is insuf- 
ficient overla p bet ween the stars in Woolf fc^Wallerstcin 



(|2005| [2006] ), and [Woolf et al.| ( |2009p and those from 



SFOCS or our CFHT samples to detect any systematic 
offsets between th e two samples. Th is is further compli 



cated by fact that Woolf et al. ( 2009 ) have very few sub- 
dwarf binaries in their sample (most of their subdwarfs 
are single stars). To test weather these low-metallicity, 
high-£ dwarfs ar e anomalous, we ob served an additional 



set of stars from Woolf et al. (2009). 



In total we observed 22 stars used in the Woolf et al.l 



(2009) calibration with SNIFS. We specifically select 
stars to cover a wide range of metallicities to get a wide 
range of C, values. We list the 22 stars in Tableland show 
them in Figure[6]as unfilled points. These 22 points form 
a clear metallicity sequence , showing that the revised £ 
can reproduce the results of Woolf et al. (20091 and that 
£ can be measured using modest resolution spectra. 

Assuming the metallicities from Woolf et al. ( 2009 ) are 
reliable and consistent with our own, we use the com- 
bined set of our binaries and the 22 additional late-type 
dwarfs from Woolf et al. ([2009) to derive the following 
relations: 



[Fe/H] = 1.55<- 
[M/H] = 1.29C- 



1.62 
1.35 



(26) 
(27) 



which we show in Figure [6] as dashed lines. The result- 
ing fits yield R^ p values of 0.58 and 0.61, respectively, 
both well above the R^ and (0.17 and 0.14, respectively). 
The RMSE values are 0.28 and 0.23, although it is no- 
tably higher for dwarfs with £ > 0.825 and lower for 
dwarfs with £ < 0.825. If we remove the 11 stars with 
[Fe/H]< -0.5 but C > 0.825, C follows a clear trend 
over the full range of metallicities covered. Interestingly, 
these 11 dwarfs appear to follow a completely different 
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Figure 6. Metallicity of the primary star ([Fc/H 
[M/H] on bottom) vs. the f para meter as defined by 
l |2007| l and using the calibration of |Lepine et al.| | |2012[ l. It nas been 
shown that the effectiveness of £ varies as a lunction of spectral 
type, so we break up our sample into four s pect ral type ranges in 
the plot. The best-fit lines from Equations |25| and |25| are shown 
as dotted lines. We add in a sample of 22 stars from Woolf et al 



(2009 open symbols) to our own wide binary sample (filled points 
(Jur measur ements of C do acc urately identify metal-poor stars 
identified in Wo olf et al.| ( [2o09^ , however, ( only identifies one of 
our wide binary stars as a subdi 



iwarf. 



sequence in [Fe/H] (or [M/H]) versus and are well sep- 
arated from their single-star, metal-poor counterparts, 
suggesting that they are unique in some way. How- 
ever, further inspection of these 11 pairs does not re- 
veal anything that could explain their discrepancy: none 
exhibit H-a significant emission (likely inactive), they 
cover a wide range of spectral types (K7-M5), and they 
have metallicities from 5 different sources (including from 
SPOCS and our own CFHT spectra). We revisit the is- 
sue of these stars in Section [8] 

7.2. J — K Metallicities 



Johnson et al. ( 2012 1 find a relation between the J — K 
and V — K colors and the me tallicity of M d warfs , based 



in part on relat ions noted by Leggett ( 1992 1 and Lepine 
fc Shara| p005| . They find a best tit reTaTIon of: 

[Fe/H] = -0.050 + 3.520A( J - K), (28) 

where A( J — K) is defined as: 



A(J-K) 



( J - K) - 0.835 

{J - K) -Y,i= <V ■ 



Ky 



V -K < 5.5 

V -K > 5.5 



and | a| = {1.637,-0.2910,0.02557}. Johnson et al 



(20121 note that this metallicity relation is only valid 
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for stars with -0.1 < A( J - K) < 0.1 and V - K > 3.8, 
but that this technique yields metallicities accurate to 
±0.15 dex. When we apply these two restrictions to our 
sample, we have 118 M dwarfs with known metallici- 
ties for their primary stars. This includes stars without 
SNIFS/IRTF spectra that were therefore not included 
earlier analyses. We find a higher RMSE of 0.20 dex, 
and an R^ p of 0.20. One possible issue is the quality 
of V magnitudes in our sample, which come from a va- 
riety of sources. However, when we remove stars with 

V — K > 5.5 unle ss they have mor e reliable V magni- 
tudes from Tycho (H0g et al. 20001, the quality of the 
fit does not change m any significant way (for stars with 

V — K < 5.5, A( J — K) is independent of V so these are 
not removed). 

We attempt to improve on the calibration and derive 
a relation for [M/H] and find: 



[Fe/H] = -0.11 
[M/H] = -0.09 



3. 14A( J -iT) 
2.UA(J-K) 



(29) 
(30) 



which results in a slightly improved RMSE = 0.19 and 
0.16, and improved R^ p = 0.30 and 0.25 for [Fe/H] and 

[M/H] respectively. These R^ p values are significantly 
larger than the R;? and values (0.08 and .09). The ma- 
j or dif ference between our fit and that of | Johnson et al.| 
(2012) is that they fix the constant term to -0.05 in order 
foTceep [Fe/H] = -0.05 at A( J - K) = 0, consistent with 
a volu me limited sample of stars from Johnson & Apps 
( |2009| ), whereas we make no such restrictions. 

7.3. K-band Metallicities 

R10 have shown that one could derive M dwarf metal- 
licities from if-band spectra using the Na I and Ca I 
lines (at 2.21/um and 2.26/im). T12 refine the calibration 
of R12 using SpeX data, and find that metallicities de- 
rived this way are accurate to ±0.12 dex. However, both 
R12 and T12 use relatively few wide binary pairs (18 and 
22, respectively), and there is significant overlap in their 
two samples. Our sample has overlap with theirs, but is 
large enough to serve as a robust check on their calibra- 
tions. Because the work of T12 was optimized for SpeX, 
we perform our test on their calibration. We follow their 
method as closely as possible (including altering our con- 
tinuum fitting procedure to match theirs). 

We find that following the calibration of T12 yields 
RMSE = 0.16 and R 2 ap = 0.69. We improve this calibra- 
tion, and find a best fit of the form: 

[Fe/H] =0.19 x EWat„ + 0.074 x EW C „ 

+2.13 x (H 2 0-K) - 3.18 (31) 

This new form yields metallicities accurate to RMSE 
= 0.14 and R ap = 0.76. The new calibration notice- 
ably improves the fit for stars with —0.3 <[Fe/H]< +0.0, 
however, both calibrations do a relatively poor job fit- 
ting the most metal-poor stars ([Fe/H] < —0.5) in the 
sample. Adding square terms improves the fit only neg- 
ligibly (AR ap < 0.02) and does not significantly improve 
the results for the most metal-poor stars. This, com- 
bined with our results from Section [6] suggests that fit- 
ting metal-poor stars requires a different set of lines, 
rather than simply higher order terms. Improvements 
may also be possible by deriving a separate calibration 



for [Fe/H]< —0.5, however, our sample has only 12 stars 
in this range, which is insufficient to derive a reliable 
calibration. 

We also find a calibration for determining [M/H] of the 
form: 



[M/H] =0.16 x EWjv a +0.039 x EW Co 
+2.29 x (H 2 0-K) - 3.04 

which gives RMSE = 0.13 and RL = 0.75. 



(32) 



0.13 and Rf p 

7.4. H-band Metallicities 

In addition to refining the calibration of RA10, T12 
derive metallicities from TJ-band spectra. The technique 
relies on the Ca and K lines in the -ff-band and a sep- 
arate H2O band defined for the ii-band (H2O-H). As 
we did with if-band metallicities in Section |7.3| we use 
our sample to test the quality of the T12 technique. As 
before, we follow their prescription, including using the 
same continuum regions to fit the continuum to a 4th 
order Legendre polynomial. 

We find the T12 calibration gives RMSE = 0.16, and 
R ap = 0.71. As before, we improve this calibration, and 
find a best fit of the form: 

[Fe/H] =0.55 x EWj? + 0.32 x EW Ca 

+1.1 x (H 2 0-H) - 2.09, (33) 

The new calibration gives an almost negligible improve- 
ment over T12; yielding RMSE= 0.14, and R* p = 0.74. 
As with the if-band metallicities, adding square terms 
improves the fit negligibly (increase in R^ p < 0.01), again 
suggesting that more lines are needed to fit the metal- 
poor stars. 

Fitting these features to [M/H] we find a best fit of the 
form: 

[M/H] =0.41 xEWjfi 0.24 x EW Ca 

+1.04 x (H 2 0-H) - 1.77, (34) 

which gives RMSE= 0.12 and R^ p = 0.71. 

8. SUMMARY AND DISCUSSION 

We present our sample of 112 late-K and M dwarfs 
in wide binary systems which we use to locate the 
most metal-sensitive features and recalibrate existing 
methods to determine late K and M dwarf metallicities. 
We combine published metallicities of 62 of the primary 
stars with 50 from our own CFHT spectra. We use 
moderate-resolution visible and NIR spectra of the late 
K and M dwarfs to identify the largest possible set of 
metal-dependent spectral features in late K to mid M 
dwarfs for each of the JHK and visible wavelength 
bands. We utilize the metallicities of the primaries to 
calibrate these metal-dependent features and obtain 
optimal relationships to estimate metallicity in M 
dwarfs. Our sample covers a wide range of spectral 
types (from K5 to M6) and 1.5 dex in metallicity. This 
enables us to search for dependencies on spectral type, 
which was previously impossible with the relatively 
small samples used. We draw 5 important conclusions 
from our analysis: 

1. It is possible to determine accurate (RMSE 
< 0.1 dex) metallicities for late-K to mid-M dwarfs 
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using modest resolution spectra (1000 < R < 2000) from 
a variety of different wavelengths. Although features 
in the K-band perform best, metallicities can be esti- 
mated from spectra of any of the four wavelength regions. 

2. Determining reliable metallicities at visible wave- 
lengths requires different calibrations depending on the 
spectral type of the star. The results are most accurate 
for K5.5-M2 dwarfs, most likely because the atomic 
lines we use are less contaminated by molecular lines 
and the pseudo-continuum is easier to estimate for these 
dwarfs. It is not known if our calibrations are applicable 
for stars later than M6. This will be the subject of 
a future investigation on metallicities for late M and 
brown dwarfs. 

3. Existing methods to determine metallicities us- 
ing H- and -RT-band spectra (e.g. those from T12) 
work well for stars of near solar-metallicity, but have 
difficulties with most metal-poor stars in our sample 
([Fc/H]< —0.5), even after applying our re-calibrations. 
Instead, determining metallicities for these stars requires 
the use of additional lines/features to improve the fit. 
This is most likely due to differences in [a/Fe], which 
are not being accurately captured by the K, Ca, and Na 
lines (Na and K are not a elements) used by R10 and 
T12. 

4. We are approaching the limits of what is possible 
with moderate resolution spectra. There is a diminishing 
return on adding additional lines to a given fit after 3-4 
features, even if there are many more metal-sensitive 
features present in a given wavelength range. Thus going 
to higher S/N or adding more wide binaries of near 
solar-metallically is unlikely to improve the calibration. 
However, improvements could probably be made by 
including later spectral types (later than M5), getting 
more [M/H] values, a larger number of more metal-poor 
stars ([Fe/H]< —1.0), or obtaining spectra with higher 
resolution. 

5. Although the Q parameter, commonly used to 
place stars into metallic ity classes, correctl y identifies 
metal-poor stars used in Woolf et al. ( 2009 ) , classifica- 
tions based on £ incorrectly identify 12 ot the 13 K/M 
companions with [Fc/H]< —0.5 as near solar-metallicity. 
This suggests that the ( parameter is sensitive to stellar 
characteristics other than temperature and metallicity 
(e.g., activity, gravity, etc.), and may incorrectly identify 
some metal-poor stars as having near-solar abundances. 

Our calibrations may be useful for both existing and 
future catalogs of M dwarf spectra. In particular, our 
calibration for visible wavelength spectra can be used 
on ex i sting catalogs of local M dwarfs such as |Lepine| 
et al. (20121 to better probe the metallicity distribution 
ot the local neighborhood especially since this sample is 
mostly early M-dwarfs, where our calibration performs 
best. Sloan Digital Sky Survey also has ~ 70, 000 vis- 



ible wavelength spect ra of M dwarfs (West et al. 2011 
Bochanski et al. 2011 ) with similar resolution to our own, 



existence of an 'M dwarf problem' ( Woolf fe West[ 2012), 
but this depends on less metallicities derived from £ pa- 
rameter. 

Although our fits have better RMSE values for [M/H] 
than they do for [Fe/H], this does not necessarily mean 
those fits are superior. In fact, the R^ p values for fits 
to [M/H] are all inferior to those calibrations done on 
[Fc/H]. The reason is that the distribution of values for 
[Fc/H] is not the same as it is for [M/H]. Some of the 
most metal-poor stars do not have [M/H] values, and 
those that do generally have higher [M/H] values due to 
large differences in a abundance (as determined for the 
primary star). 

We confirm the claim of Johnson et al. ( 2012[ ), that one 
can get approximate M dwarf metallicities using J — K 
versus V—K colors. However, the technique has a limited 
range of metallicities (—0.4 <[Fc/H]< +0.2) and is only 
accurate to ~ 0.2 dex. Thus this technique is probably 
best used in special applications, such as biasing a planet- 
search towards metal-rich M dwarfs. 

One possible explanation for the poor performance o f £ 
on our sample compared to that of Woolf et al. ( 2009 ) is 
the presence of unresolved binaries. It is likely that most 
wide binaries form as higher order systems (Kouwen 



hoven et al. 2010). Thus many of our wide pairs may 
include unresolved M+M dwarf pairs. There is evidence 
of radius inflation in low-mass eclipsing binary systems 



( |L6pez-Morales|2007||Irwin et al |2011 |Kraus et al.|2011 I 

and may also be cooler th an their single star counter- 
parts (Boyajian et al. 2012). Further, atmospheric mod- 
els indicate that the Ti05 and CaH2/CaH3 indices, on 
whi ch C is based, are sensitive to tempe rature and grav- 
ity ( Jao et aL| |2008 Allard et al. 2011). However, none 
of our most metal-poor companions show H-a emission, 
whereas radius inflation in tight binaries is usually asso 



dated with high chromospheric activity jLo pez-Morales 

T 



2007 Kraus et al. 2011 



Stassun et al. 20 



2l. Additional 



which could be used in conjunction with our calibrations 
to map out the metallicity distribution of the sample. 
Work has already been done in this area to confirm the 



metallicities ot M-dwart with known multiplicity (e.g., 
low-mass eclipsing binaries and spectroscopic binaries) 
are needed to confirm if this is the source of the discrep- 
ancy. 

Another complication is the possibility of having false 
binaries (chance alignments) in our sample. We can 
estimate the number of in terlopers by cross referenc - 
ing our sample with that of |Tokovinin fe Lepine| ( |2012[ ). 
Tokovinin & Lepine (2012) calculate the probability that 
stars with, commiserate proper motions are actually phys- 
ically associated wi t h eac h other (P p h ys )- Although 
Tokovinin & Lepine (2012) caution that their probabili- 
ties are purely based on models (and therefore only ap- 
proximate), the numbers can be used to give a rough 
estimate of contamination from chance alignments. By 
summi ng up P„h V a values fo r all of our binaries included 
in the Tokovinin & Lepine ( 2012[ ) sample we find that 
> 90% of our binaries are physically associated with each 
other. However, some of the pairs with low -Pphys values 
have parallax information for both the primary and com- 
panion that are consistent with each other. If we assume 
pairs with consistent parallaxes have P p hy S — 1 and re- 
peat our calculation, we find that 94% of our binaries 
are physically associated with each other. Although our 
metal-poor stars ([Fe/H]j-0.5) tend to be more distant, 
and therefore only 3 of the 13 are listed in Tokovinin fc 
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Lepine (2012), two of them have -Pphys > 93% (the other 
has -Pphys = 72%). Three more of our [Fe/H]< —0.5 stars 
have parallaxes that are consistent with the primary to 
la, indicating that even our metal-poor stars are almost 
all physical pairs. 

We do not claim to have identified every single metal- 
sensitive feature at the resolution of our spectra, how- 
ever, the nature of our analysis means that it is unlikely 
that we missed any of the most useful ones. We per- 
form a rough test on our recovery rate by introducing 
artificial metal-sensitive lines of various usefulness and 
repeating our analysis. Specifically, we select a sample 
of the most metal-sensitive features (those used in Equa- 
tions (8)-(17)) and insert them elsewhere in the spec- 
trum of the stars. When moving features from NIR to 
visible wavelengths we convolve the lines with a Gaus- 
sian profile to reduce the resolution of the feature (fea- 
tures moved into the NIR are not changed). We then 
repeat our analysis as described in Section [5] We find 
that metal-sensitive features arc sometimes not identi- 
fied when they are placed on very strong telluric fea- 
tures, blueward of 0.4/im (where the S/N is very low), 
or when they overlap with other strong features (e.g., 
the Mg lb line) that make clean measurements difficult. 
We also note that features identified as metal-sensitive in 
the NIR appear less metal-sensitive (although they are 
still identified) when placed in visible wavelengths; this is 
likely due to the lower resolution and/or difficulties mea- 
suring features that are convolved with strong molecu- 
lar lines in the visible. In spite of these exceptions, we 
still recover > 88% of the lines on average, and > 93% 
when we exclude telluric regions and low S/N regions of 
the spectrum. This indicates that our analysis is quite 
robust, and that expanding on our findings will require 
observations later-type stars (past M5), more metal-poor 
stars, or higher resolution visible spectra. 
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